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CN [ Abstract 

Given a convex set and an interior point close to the boundary, we prove the existence 
of a supporting hyperplane whose distance to the point is controlled, in a dimensionally 
\Q , quantified way, by the thickness of the convex set in the orthogonal direction. This result 

has important applications in the regularity theory for Monge- Ampere type equations arising 
in optimal transportation. 

< 



1 Introduction 

In this note we establish an estimate which quantifies the dimensional dependence of the claim 

that corresponding to any (interior) point near the boundary of a convex set, is a supporting 

hyperplane much closer than the thickness of the set in the orthogonal direction. The main 

\q • interest of our estimate (see Theorem II. ip is that it allows us [1] to extend — for the first time 

- a Holder continuity result of Caffarelli [21I3J concerning optimal transportation of bounded 

measurable densities from the quadratic cost function of Brenier [TJ, to the more general cost 

functions considered by Trudinger and Wang [7]. Caffarelli's regularity result has well-known 

connections to convex geometry (see [5]), and the importance of its extension to more general 

7— I \ optimal transport problem is highlighted in [8]. 

Our theorem below is purely geometric, elementary to state, and non-trivial to prove; it 
may well be of independent interest. To emphasize this possibility, and the fact that it does 
K> , not rely on any auxiliary concepts arising from the intended application [3], we establish it in 
this separate manuscript. By so doing, we hope to ensure its accessibility to convex geometers 
who may have no interest in optimal transportation, as well as to its primary target audience, 
which consists of researchers interested in the regularity of optimal mappings (or equivalently, 
of degenerate elliptic solutions to the associated Monge- Ampere type equations). 

Let us start by recalling that a convex body S in the n-dimensional Euclidean space R n refers 
to a compact convex set with non-empty interior. A well-known result of Fritz John [6], often 
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called John's Lemma, shows every convex body can be translated so that it contains an ellipsoid 
E centered at the origin, whose dilation by factor n contains the translated copy S of S: 

E C S CnE. (1.1) 

The constant n is shown to be sharp by the standard simplex. After this translation, i.e. when 
(jl.ip holds, we call S well-centered. We restrict our discussion to only well-centered convex 
bodies, but this does not cost any generality. 

For any point y € S near the boundary of a well-centered convex body, we claim it is 
possible to find a direction in which the boundary of S is much closer than the thickness of S 
in the same direction. More precisely, we claim it is possible to find a line L through the origin 
whose intersection with S is large relative to the distance of y to a hyperplane outside of S and 
orthogonal to L. Here orthogonal refers to the ambient Euclidean inner product, so that the two 
distances being compared are measured along line segments parallel to L. The following theorem 
quantifies the dependence of their ratio on the proximity of y to the boundary, as reflected in 
the degenerating factor s 1 ' 2 " in (jl.2p below. When y approaches the boundary of S, the ratio 
of the two distances becomes more and more exaggerated, algebraically fast with respect to the 
separation of y from the boundary (but whose algebraic power decays exponentially fast in high 
dimensions). Here s measures the separation of y from the boundary in the Minkowski gauge 
of S — which of course is equivalent to any other norm on R n . In what follows however, 'dist', 
'diam' (and orthogonality) always refer to distance and diameter with respect to the Euclidean 
norm. The application [1] requires only the special case sq = l/(2n). 

Theorem 1.1. Let S C R n be a well-centered convex body, meaning (jl.ip holds for some 
ellipsoid E centered at the origin. Fix < s < so < 1. For each y £ (1 — s)dS there exists an 
hyperplane P supporting S such that 

dist(y, P) < e(n, s )s 1/2n ~ 1 diam(P ± n 5). (1.2) 

Here P denotes the (unique) line orthogonal to P passing through the origin, and c(n, so) is a 
constant depending only on n and sq, namely c(n, sq) = r?' 2 {n — |) I - — (L /2 n I 

Remark 1.2. For n = 1, the constant c(1,sq) is sharp and (|1.2p becomes an obvious equality 
in that case; we have not investigated sharpness of c(n, So) or of the power 2 1_n in higher 
dimensions. 

The key point of the estimate (jl.2p is that the ratio dist(y, P)/diam(P J - n S) goes to as 
s — > in a "uniform way" , independent of the shape of S. Observe that if s = we can choose 
P to support S at y € dS, but for s > it is less obvious how to choose P (and hence L). The 
difficulty for proving this estimate is on the arbitrariness of the convex body S. For example, 
if S is the round ball, then the estimate (jl.2p becomes trivial (indeed, one may even replace 
c(n, sq)s 1 ' 2 with s, and the supporting hyperplane P shall be chosen to be orthogonal to the 
vector y). For a general convex bodiy S, there are three natural ways to try to generalize such 
choice: (i) to choose P orthogonal to y, (ii) to choose P supporting S at the intersection of the 
half line Oy with dS, or (iii) to choose P closest to y. However, in all these three cases it is not 
hard to find counterexamples (some family of degenerating thin convex bodies) showing such 
rather natural choices of hyperplanes do not work, namely, not yielding a uniform convergence 
to of the ratio dist(y, P)/ diam(P J - n S) as s — > 0. To prove Theorem II .11 we find an algorithm 



which allows to choose appropriate P by an inductive recursion, reducing the dimension of the 
task confronted at each step. 

One of the reasons why the estimate (|1,2|) is nontrivial is that the set of convex bodies is 
not compact. A common and powerful way to deal with such non-compactness is to use John's 
Lemma |6] (see (jl.ip ). to renormalize the convex bodies via afhne maps so that the resulting 
shapes become roughly close to the round ball (with a uniformly bounded scale factor). This 
way, one can easily derive some estimates for quantities that are affine invariant. One such 
example is the classical Alexandrov estimates for the Monge-Ampere measure associated to 
sections of convex functions: see, for example [5]. However, in our case the inequality (11.2p 
involves orthogonality with respect to the fixed Euclidean norm, which is not affine invariant. 
Therefore, we cannot derive the estimate (jl-2j) by applying John's Lemma. 

The remainder of this paper is devoted to the proof of Theorem ll.lt which is completely 
elementary though quite nontrivial. 

2 Proof of Theorem 11.11 

Let S C R n be a well-centered convex body. Let the positive numbers a ,■■■ ,a n € -R+ 
denote the lengths of the principal semi-axes of the inner ellipsoid E of Fritz John (jl.ip . One 
can regard these a J 's as the coordinate components of the vector a = [a , • • • ,a ra ). In the 
following, superscripts will be used to denote such coordinate components for vectors, and for 
all the other cases they will mean powers (we believe this should not create either confusion or 
ambiguity) . 

Use these principal axes to choose coordinates, with the origin at the center of E. We still 
have the freedom to choose the order in which these axes are enumerated, which we shall exploit 
especially at (12.9p . In these coordinates E is represented as 



c n % o ^ 

E= , = (, 1 ,..,/) e ^|£g) <1 . 



The rectangle 

R n = {x = (x 1 ,- ■ ■ ,x n ) € R n | \x l \ < na\ i = 1, • • ■ ,n} 

circumscribed around the outer ellipse nE will also play a crucial role. Observe -^R n C S C 
R n . In particular, S is comparable (in size and shape) to R n . 

2.1 Initial step in the recursive algorithm 

Fix so € (0, 1), and given y £ (1 — s)dS, < s < sq, let p n be the intersection of the half line 
Oy with dS. We pick a tangent hyperplane P n (which may not be unique) to S at p n . Using 
similar triangles, we deduce that 

dist(y,P n ) 



dist(0,P n ) 
Two alternatives 

We consider the following two exclusive cases: 



Favorable case: Suppose we are lucky enough that 

>s 1 ' 2 . (2.1) 



diam(P^ n R n) _ a/2 



2dist(0,P n ) 

Then, the choice P = P n leads to the desired result (jl.2p since diam(P r ^ni? n ) < n 3 ' 2 diam(P^-n 
S), thus 

dist (y,Pn) 3/2 dist(y,P n ) dist(0,P n ) 



diam(P^ n S) ~ dist(0, P n ) diam(P^- n R n ) 

<n 3 / 2 - '' 



7,3/2 

n 1/2 



2aV2 

s 1 



Unfavorable case: If the convex body S is very thin, or equivalently if the outer rectangle 
Rn is very thin, then the ratio diam(P^ L n R n )/2dist(0,P n ) can be much smaller than s 1 ' 2 , in 
which case (|2,ip fails. For such situations, we now describe a recursive algorithm which shows 
that whenever (|2.ip is violated, after at most (n — l)-steps it is possible to find an alternative 
hyperplane P (generally different from P n ) which fulfils the desired conclusion (jl.2p . 



2.2 Notation in the recursive algorithm 

The basic idea of the following recursive algorithm is to repeat the previous two alternatives 
in the inductive steps, with decreasing dimension. Since this is a finite dimensional situation, 
such algorithm should terminate, and we show it does so yielding the desired result (jl.2p . One 
of the key points of the argument is to choose the right geometric configuration. This requires 
in particular some careful choice of the terms 7«(s), 5i(s), and cq, as we fix the notation below. 
Let i, k G {1, . . . , n) and s G [0, sq]. We define 



S i (s):=(2i-l)s 1 / 2i 



These satisfy the following key relations: 



*i+l(a) > Si(s) + 2 7i (s); 
Si(s) 



(2« - lH(s). (2.2) 

ii\ s ) 

We use coordinates (x , • • • , x +1 ) on R + . Define the projections 7T& : R + — >■ R by 

7T k (x\--- ,X k ,X k+1 ) = (x\--- ,X k ). 

Observe that each ti^ is determined by the choice of the coordinate axis for x for R + , 
and such choice will be made individually at each step of the recursive algorithm. This is an 
important point to remember throughout the proof. For k < n, define the rectangles Rk in R ' 
inductively as dilated projections of R n : 

Rk = CQTTk(Rk+l)- 



Here, the dilation factor cq (with respect to the origin 0) is given by 

7*0) 



Co := max 

0<s<so; l<i<n 



1 + 2 



1-7*00 



(2.3) 



-i±M^ (24 ) 

1 - (so) 1/i! 

where the monotone dependence of ji(s) on both s < so and i < n has been used. These 
rectangles can also be written as 

R k = { x = (a; 1 ,... ,x k ) G i? fc | \x l \ <c^~ k na\ \<i<k}. (2.5) 

Let Q k (and Q^) be the parallel hyperplanes in R which form the boundary of (and bisect) 
R k orthogonally to the x fe -axis: 

Q± = { x = ( x \ . . . , x k ) e R k | x k = ±$- k na k }; 
Q° k = {x = (x 1 ,--- ,x k ) eR k | x k = 0}. 



Set y n = y = (y 1 ,...,y n ) and define its projections recursively 

Vk '■= n k °n k+1 o ••• o-7T n _i(y) = (y 1 ,...,y k ). 

Since cq > 1 it is clear that y k G R k - 

In the following we will define some hyperplanes P k C R inductively. For such a hyperplane 
Pk C R use Pfc to denote the extension of P k to the hyperplane in R n parallel to the x k+1 - 
through x n -axes, i.e., which satisfies 

7Tfc O ••• OTT n „i(P fe ) = P k . 

2.3 The structure of the recursive algorithm 

Recall the point p n = Oy fl dS and the supporting hyperplane P n , with p n € P n , which are 
given in the initial step of the recursive algorithm (Section 12. ip . Let us describe the recursion 
in which we use (p k ,P k ) to define (p k _i, P k -i)- We first list three conditions that are required 
at each step of the recursive algorithm. For the i step (here i = n — k + 1) assume that 

A fc : p k e P k n R k , y k G [0, p k }; 
dist(y k ,P k ) 
k dist(0,P fc ) " n k+n h 

C k : P k does not intersect the interior of S. 

Notice that for k = n, either the favorable situation (|2,ip holds, in which case there is nothing 
further to prove (and so no need to proceed to the next step), or else A n ,B n ,C n are satisfied 
by our initial choice of p n and P n . 



Starting from k = n, we shall decrease k one step at a time until the algorithm terminates. 
Whether or not the recursion terminates at a given value of k is determined by the following 
dichotomy: 

Case I (favorable case): -^ fc k > J n ~k+i{s); (2.6) 

2dist(0,P fc ) 

rliam ( F*^- D /? ^ 

Case II (unfavorable case): 2d' tfO P ) < 7n-k+i(s). (2.7) 

If Case I holds for some value of k we shall discover we are in a favorable situation — analogous 
to (|2.ip — which allows us to terminate the recursion and obtain the desired result (|1.2p . On the 
other hand, if Case II holds for the given value of k, we shall see we can use (p k , P k ) satisfying 
the inductive hypotheses A& , B^ and C& , and the condition (|2.7[) . to construct (p fe _ 1 ,Pfe_i) 
satisfying Afc_i , B^_i and Cfc_i . We then decrease fc and proceed to the next step of the 
recursion. In the worst case the recursion continues until k = 1, and we find (pi,P\) satisfying 
Ai , Bi and Ci . In this case we show in the last section below that the desired result (|1.2p can 
again be obtained, to complete the proof of the theorem. 



2.4 Case I, ( 12. 6 ft holds for some k > 2: the recursion terminates with the 
desired result. 

As soon as we reach some k > 2 for which the condition (j2.6h holds, we stop the recursion. 
We now show in this case the desired result (|1.2p follows. Here, the assumptions B& and C k are 
crucial. Since, by the construction of y k and Pk, dist(y,P k ) = dist(y k , Pk) and dist(0, P&) = 
dist(0, Pfc), recalling (|2.5p we get 

dist(y,P fc ) dlst{y k ,P k ) dist(0,P fe ) 



c"~ fc diam(P fc ± n R n ) " dist(0, P k ) diam(P fe 1 n P fe ) 

< ^" fc+l( ; } (bydZ6l)andB fc ). (2.8) 

27n-fc+l( s J 

Let H k be the half-space containing 0, with dH k = P k . Notice that S C H k by assumption C k . 
Thus, translating P k toward S, one can find a hyperplane P supporting S, which is parallel to 
P k . Since dist(y, P) < dist(y, P k ) and diam(P ± nS) = diam(P^ n5)> n~ 3 / 2 diam(P fc - L n P„), 
from ([272]) and ([275]) we have 

diSt(y,P) < 3/2 n-ifc ^n~fc+l(s) 



diam(P^n5)" ° 2 7n _ fe+1 (s 



„3/2^-fc (n _ fc+ l )a l/2»-«H-l 



1TJ-1 



(recall that k > 2), which gives the desired result (|1.2p . 

2.5 Case II, (12.71) holds for 2 < k < n: the recursion continues. 

A remark before we proceed: in the following argument, we assume that any claimed intersections 
between affine subspaces such as lines and (hyper-)planes actually exist and have the expected 
(i.e. generic) dimension. This costs no generality for two reasons: 



• To avoid parallelism we can perturb if necessary (i.e. rotate and/or translate slightly) the 
afhne subspaces. 

• We will obtain estimates which are not sensitive to small perturbations, so the estimates 
also hold without the perturbation. Moreover, the obtained bounds then imply that the 
claimed intersections do indeed exist. 

Similarly, we can also avoid, if necessary, the cases where some lengths and/or distances degen- 
erate to zero. 

Some preliminaries and the construction of Pk-i 

To define P k —x and Pk-\ we set-up preliminaries. Let us first consider the point r^ = 
( r fc' ' ' ' ' r k) e -^ defined as the closest point on Pk to the origin, so that dist(0, r^) = dist(0, Pk). 
Note r k is outside the rectangle Rk, since otherwise, 

diam(P fc ± n R k ) > 
2dist(0,P fc ) " 

contradicting assumption (|2.7p . Let r^ = ((r^) 1 , • • • , (r~^) k ) = [0, r^] C\dRk denote the intersec- 
tion point of the ray through r k with whichever of the 2fc faces of this rectangle it intersects. 
Without loss of generality, suppose the axes are enumerated so that the intersection occurs on 
the face of Rk contained in Q^ . (Observe that this choice of coordinates affects the definition 
of 7Tfc_i.) Then, 

(4t = Co~ fe a fc > (2.9) 

holds. Because dist(0,r^~) = ^ diam(P^- n Rk) and dist(0,T"fc) = dist(0,Pfc), we have 

dist(0,r+) = diam(P fc ± n R k ) 

dist(0,r fc ) 2dist(0,P fc ) ' l ' ' 

Now, define the hyperplane Pk-i in R ~ l by 

Pk-i-irk-iiPk^Qk). (2.11) 

Verification of Cfc_i 

Before proceeding further, let us verify that Cfc_i follows from C^ as a consequence. 

Since r\ = ^/kfcl 2 ~~ l 7r fc-i( r fe)| 2 ) by a simple geometric argument (see Figure 1) the con- 
struction above yields T\.-\ = ^k-i( r k) (recall that r k -\ is defined as the closest point on Pk-i 
to the origin), with A = (|r fc | 2 + c^c^rD/^k-iirk)? > 1. 

For each x = {x , . . . , x n ) =: x n from the interior of S, let Xk '■= ftk ° ^fe+i • • • ° vr n _i(a;) = 
(x 1 ,...,x k ). To verify Cj..! we need to show that Xk-\ ■ r k -i < |f/c-i| 2 - Without loss of 
generality assume that Xk-\ ■ T^—x > 0. Hence, since A > 1, 

Xk-i ■ r k -i < X(x k ■ r k + a k r k ) 
< \(\r k \ 2 + a k r k k ) 
< kfc-il 2 , 

where the first inequality follows from x > — a , the second from C k ■, and the third from cq > 1 
and the definition of A. This yields C^_i as desired. 



c?-* a k 



r n-k n k 
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K-ll 



Figure 1: By a simple argument based on similar triangles, one can easily compute |7Tfe_i(rfc)| 



in terms of IrJ, d\ 



~ k a k , 



k- 



and |7r fc _i(r fc ^ 



Construction of p k _ 1 



yk-l 



We will now find p k _\ € R so that A^_i and Bfc_i are satisfied. To define p k -i, consider 
the two-dimensional plane T k C R , generated by the x fc -axis and the half-line 0p k (which is 
the same as 0y k ). (We perturb p k slightly if necessary to ensure it does not lie on the x -axis.) 
Since dim.T k + dimP^ = k + 1, the affine intersection L k := T k n P k C -R contains at least a 
line; it contains at most a line since S T k \P k . Notice that the line L k passes through the 
point p k and the hyperplane Q7, as in Figure 2. 

We define 



Pk-i '■= Kk-i(L k nQ k 



Notice that 



(2.12) 



(2.13) 



In particular, to verify A^_i we only need to check p k _\ € Rk-i- 

Preparation before verifying Afc_i and Hk—i 

To verify A k -\ and Bfc_x we first find a few relevant points on the plane T k . What we are 
going to describe is summarized in Figure 2. 

The line L k intersects with the three parallel hyperplanes Q k ,Q k C R ■ Denote the three 
intersection points by 



Ik 

4 



L k nQ+; 
Lk nQ fc ; 
L k r\Q^. 



Notice that p k € [q%, q k ] and 7r fe _i(q fe ) = p k _ v 

Let l k denote the closest point on L k to the origin 0, and let t k G T k be the orthogonal 
projection of r k on the plane T k . Notice that, since r k is the closest point in P k to and 
L k C P k , the orthogonal projection of the ray 0r k to T k is the ray 0l k and \r k \ < \l k \ (to see 



^-i=%i(e:np,) 




Ql n p, 



Figure 2: The geometric argument behind the construction of p k _ l . Observe that, for s small, 
this figure (as well as the other figures) is not very "realistic", as R k should be very thin in 
the "horizontal" directions and Pk-i very close to Rk- However, even if the proportions are not 
respected, this figure should help the reader to follow the argument described below. Note that 
one can construct convex sets for which y k and p k may lie below Q k . However, although the 
picture is just indicative, our proof is purely analytic and works independently of the position 
of y k and p k with respect to Q k . 
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Figure 3: Some useful similitudes. Note that, analogously to what observed in the comment to 
Figure 2, it maybe also possible that p k (resp, 7Tfe_i(p fe )) be located between q° k and q k (resp, 

Pk-l)- 

this, one may consider the plane passing through 0, r k , and Ik, and observe that it cuts Lj- 
orthogonally) . In particular, the point t k belongs to the line segment [(),/&]. Moreover, since 
Tfc contains the x fc -axis, t\ = r%. Hence, by our assumption (|2.9p . £& belongs to the region over 
Q^, namely t\ > Cq~ a k , as (therefore) does Ik- Consider the point 

[<M fc ]nQ+ 



it- 



with (li 



(rtY 



ji—k„k 



a . Since (as we observed above) t k = rjj, the triangles A(£jfc,0,r* 



and A(Z7~,0,r7~) are similar. Thus 



dist(0,Z+) ^ dist(0,Z+) _ dist(0,r+) 



< 



dist(0,Z fc ) ~ dist(0,t fc ) dist(0,r fc ) 
Here, the last inequality follows from (|2.10j) and (|2.7p . 



(2.14) 



Now the triangle A(0, Ik, qy f° rme d by the three points 0, l k and q® has a right angle at Zj 



fe-i 



This entire triangle projects to a line segment in R , with the projection 7Tfc_i preserving the 
order of points and ratios of distances along the edges of A(0, Ik, q\) as in Figure 3 — a fact we 
shall continue to use subsequently. 

Similarity to A(Z7~, Z&, qt) combines with (|2.14p to yield 



dist(7T fc _i(gg),7r fc _i( 9 +)) _ dist(q°,g+) _ dist(0,Z+ 



dist(7r fc _ 1 (q°),vr fc _i(Z A .)) dist(q°,Z fc ) 
Now we are ready to verify Afc_i and ~Bk-i ■ 

Verification of A^_! 

Claim 2.1. Definitions ([2TT]) - (f2T2"]) imp/y A fc _i . 



dist(0, Zj 



< 7 n _ fc+ i(s). 



(2.15) 



Proof. By (J2.13J) it is enough to show that p k -\ € Rk-i- This should be clear from the geometric 
construction. Here the factor cq in the definition Rk-i = Co^k-i(R-k) plays a crucial role. To 
give details, first note that Afe implies 

7Tfe-l(Pfe) G 7Tfc_l(jRfc). 



10 



From p k G [q£, q k ] and p k _ x = n k _ 1 {q k ), we have 

dist(0,p fe _ 1 ) = dist(0,7T fc _i(p fc )) + dist(7r fe _i(p fc ),P fe _i) 

< dist(0,7r fc _i(p fe )) + dist(7r fc _i(q+),p fe _ 1 ) 

= dist(0,vr fc _i(p fc )) + 2dist(7r fe _i(q+),vr fc _i(g^)). 

Here, to bound the last line, observe that from (|2.15j) . 

dist(7rfe_i(<z£),7Tfe_i(qf°)) < - — n ~ k+1 dist(7r fc ^i(Z fc ),7r fc „i(g+)). 

t — 1n-k+l{S) 

From the geometry of the right triangle A(0, l k , q k ) C T k , 

dist(vr fe _i(« fc ),7rfc-i(9fc)) ^ dist(0,7TA._i(p fe )). 

Recalling that Pfc_i is parallel to n k —i(p k ), combining the preceding four displayed statements 
with (|2.3h yields p k _i G coTTk-i(Rk) = -Rfc-i as desired. This completes the proof of Claim |2~T1 

□ 



Verification of Bj._i 

Claim 2.2. Definitions ([2TTT]) - ([2712]) imp/y B fc _i . 

Proof. Since 2/fc_i G [0,p fc „ 1 ], from similarity 

dist(y fc _ 1 ,P fc _i) dist(y fc _ 1 ,p fc _ 1 ] 



dist(0,i\_i) distCO.Pfc.!) 

To bound the latter: 

dist(Vfc-i,J»fc-i) 

= dist(7r fe _i(y fc ),7r fc _i(p fc )) +dist(vr fc _i(p fc ),p fc _ 1 ) 

< <5 n _ fe+ i(s) dist(0,vr fc _i(p fc )) + dist(vr fc _i(p fc ),p fc _ 1 ) (by B fc and similarity) 

< 5 n _ fe+ i(s)dist(0,p fc _ 1 ) + dist(vr fe _ 1 (p A .),p fc „ 1 ) (by [0,n k -i(p k )] C [0,p fc _x]) 

< <5 n _ fc+ i(s)dist(0,p fc _ 1 ) + dist(vr fc _i(q+),p fc _ 1 ) (by [7T fe _i (p^p^] C brfc-ite^Pfc-i]) 

< 5 n _ fc+1 (s)dist(0,pj : _ 1 ) + 27 n _ fc+1 (s)dist(7r fc _ 1 (Z fe ),7r jfe _ 1 (g^)) (by ([US])) 

< (<J„_ fc+ i(s) + 2 7n _ fc+1 (s))dist(0,p fc „ 1 ). (by [7r fe _i(Z fe ),7r fe _i(g^)] C [0,p fc _!]) 



< ^n-fc+2(s) 



Therefore, by 02721) . 

distd/^LPfc-i) 
dist(0,P fc _!) 

which is the desired result. This completes the proof of Claim l2~2l □ 

We have shown that the assumptions Afc_i , T5 k -i and Cfc_i are satisfied, therefore we can 
continue the recursion until we arrive at Case I, where we get the desired result (jl.2j) . or else, 
at worst, arrive at the following scenario. 
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2.6 Final remaining possibility: the recursion reaches k = 1. 

Suppose that this recursive procedure does not stop before we find Pi,Pi € R 1 satisfying 
Ai , Bi and Ci and decrease k from 2 to 1. We now show the desired result can be established 
in this case. Writing p 1 = p\ to emphasize that we are now dealing with 1-tuples, Ai yields 
P 1 = { Pl } c Ri = {x e R I \x\ < do' 1 a 1 }, so 2dist(0, Pi) < diam(P 1 ± n Ri). Therefore, 

dist^^Pi) dist(y 1; Pi) 



diam(P^ n Rt) ~ 2dist(0,Pi) 

<8 n (s)/2 (byBi). 

Since C\ guarantees Pi is disjoint from the interior of S, we can argue exactly as Case I to 
show the supporting hyperplane P of 5 parallel to Pi satisfies 

dist(y,P) / 3/2. dist(y,Pi) 



< n° 



diam(P ± n 5) " diam(P 1 J - n R n ) 

,3/aj-i dist(yi,Pi) 



n 



diam(P 1 - L nPi) 
C n*' 2 (%- x 8 n (8)l2 

= n^ 2 c n ^ 1 {n-\)s 1 l 2n -\ 
the desired result (jl,2p . 

References 

[1] Y. Brenier, Polar factorization and monotone rearrangement of vector-valued func- 
tions. Comm. Pure Appl. Math., 44 (1991), 375-417. 

[2] L.A. Caffarelli, Some regularity properties of solutions of Monge Ampere equation. 
Comm. Pure Appl. Math., 44 (1991), no. 8-9, 965-969. 

[3] L.A. Caffarelli, The regularity of mapping with a convex potential. J. Amer. Math. Soc, 
5 (1992), 99-104. 

[4] A. Figalli, Y.-H. Kim and R. J. McCann, Holder continuity and injectivity of optimal 
maps. Preprint. 

[5] C. Gutierrez, The Monge-Ampere Equation. Birkhauser, 2001. 

[6] F. John, Extremum problems with inequalities as subsidiary conditions. Studies and 
Essays Presented to R. Courant on his 60th Birthday, January 8, 1948, 187-204. Inter- 
science Publishers, Inc., New York, N. Y., 1948. 

[7] N. Trudinger and X.-J. Wang, On the second boundary value problem for Monge- 
Ampere type equations and optimal transportation. Ann. Sc. Norm. Super. Pisa CI. 
Sci. (5) 8 (2009), no. 1, 143-174. 

[8] C. Villani, Optimal Transport, Old and New, Grundlehren des mathematischen Wis- 
senschaften [Fundamental Principles of Mathematical Sciences], Vol. 338, Springer- 
Verlag, Berlin-New York, 2009. 



12 



